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Abstract: The partition function on the three-sphere of N = 3 Chern-Simons-matter theories 
can be formulated in terms of an ideal Fermi gas. In this paper we show that, in theories with 
N = 2 supersymmetry, the partition function corresponds to a gas of interacting fermions in one 
dimension. The large N limit is the thermodynamic limit of the gas and it can be analyzed with 
the Hartree and Thomas-Fermi approximations, which lead to the known large N solutions of 
these models. We use this interacting fermion picture to analyze in detail N = 2 theories with 
one single node. In the case of theories with no long-range forces we incorporate exchange effects 
and argue that the partition function is given by an Airy function, as in M = 3 theories. For the 
theory with g adjoint superfields and long-range forces, the Thomas-Fermi approximation leads 
to an integral equation which determines the large N, strongly coupled i?-charge. 
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1. Introduction 

In the last years, the use of localization techniques in the study of Chern-Simons-matter (CSM) 
theories has led to a new class of matrix models which generalize the matrix models for pure 
Chern-Simons theory introduced in [22, 24]. These models compute partition functions and 
Wilson loop correlators on the three-sphere and other compact three-manifolds. They were first 
introduced in [19] for CSM theories with N > 3 supersymmetry, and the result was extended to 
theories with N = 2 supersymmetry in [14, 9], see [26] for a review and a list of references. The 
study of these matrix models has led to many interesting results. It has provided precision tests of 
the AdS4/CFT3 correspondence as well as beautiful field-theoretical results on superconformal 
field theories in three dimensions. It is therefore interesting to find efficient ways to analyze 
these models, in particular in the large N limit, where they make contact with superstrings and 
M-theory. 

The matrix model corresponding to ABJM theory was solved in the 't Hooft expansion, 
at all orders in 1/N, in [27, 6], by using techniques developed earlier in random matrix theory 
and topological string theory. However, these standard techniques seem to be of limited use for 
models with N < 3 supersymmetry (see however [5, 36] for planar solutions of some of these 
models by using those techniques). For more general theories with N = 3 supersymmetry, [10] 
developed a powerful technique to study the so-called M-theory limit, where N is large and the 
coupling constants are fixed. This technique can be generalized to M = 2 theories and it has 
made possible to understand the leading, large N limit of the free energy of these matrix models, 
see [15, 29, 4, 17, 16] for a non-exhaustive list of examples. 

In [28], a different framework was introduced to analyze N = 3 theories, based on a refor- 
mulation of the matrix model as the partition function of an ideal Fermi gas. This reformulation 
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has various virtues: first of all, it gives an elementary derivation and a nice physical picture of 
the iV 3 / 2 behavior of the free energy of these models in terms of free fermions. Second, it makes 
possible to compute all 1/N corrections to the free energy. This leads to a simple derivation of 
the Airy function result for the partition function of ABJM theory [7]. In fact, it was conjectured 
in [28] that in theories displaying the iV 3 / 2 behavior of the free energy, the partition function at 
all orders in 1 /N is given by an Airy function, and the Fermi gas reformulation makes possible to 
establish this conjecture for many N = 3 theories. The Fermi gas picture also leads to powerful 
results for the Wilson loop vevs of ABJM theory at all orders in 1/N [20], and it is natural to 
see how it can be generalized and/or used in other situations. 

In this paper we generalize some of the results of [28] to CSM theories with M = 2 su- 
persymmetry. We find that the partition function on the three-sphere of these theories can be 
reformulated, via the matrix model of [19, 14, 9], as the partition function of an interacting Fermi 
gas. For theories with M > 3 supersymmetry, the Fermi gas is non-interacting but there is a 
non-trivial external potential, and all the physics reduces to one-body physics. In the case of 
M = 2 supersymmetry, we have in general two-body and higher body interactions. Interacting 
quantum gases are notoriously hard to analyze, and this one is not an exception. However, in 
the thermodynamic limit (which corresponds to the large iV limit of the matrix model), one can 
use the Hartree/Thomas-Fermi approximation. We show that this approximation leads to the 
large N treatment of the N = 2 theories proposed in [15, 29, 4]. Unfortunately, a systematic 
understanding of the 1/N corrections to this leading, large N result seems difficult in general. In 
the cases of flavored theories with one node first considered in [3, 13], it is in principle possible 
to go beyond the Thomas-Fermi approximation, and we compute the next-to-leading correction 
due to exchange effects in the interacting Fermi gas. We also study a theory with one node but 
with long-range interactions between the fermions, namely the CS theory with adjoint multiplets 
considered in [8, 30, 33, 31]. In this case, the long-range forces lead to a qualitatively different 
physics for the interacting fermions, which can be still studied in detail at large N by using the 
Thomas-Fermi approximation. We obtain in particular an integral equation which determines 
the large N, strongly coupled i?-charge of the multiplets. 

The organization of this paper is the following. In section 2 we review the construction of 
matrix models for J\f = 2 CSM theories, we show how to interpret them in terms of interacting 
fermions, and we explain how to obtain the large N solution via the Hartree/Thomas-Fermi 
approximation, reproducing in this way the density functional approach of [10, 15] and other 
papers. In section 3 we study in detail one node, flavored theories. We analyze them in terms of 
an interacting Fermi gas, we solved them with the Thomas-Fermi approximation, recovering the 
solution of [15], and we calculate the next-to-leading, exchange correction to this leading result. 
We argue that the partition function is an Airy function. In section 4, we study from the point 
of view of the interacting Fermi gas a one-node theory with long-range forces, namely CS theory 
with adjoint multiplets, and we derive an integral equation for the large N limit of the R-chage. 
This allows us to re-derive, and partially improve, some results in [30] on this theory. Finally, 
we conclude with some remarks and open problems. 

2. J\f = 2 Chern— Simons— matter theories as interacting Fermi gases 
2.1 Matrix models for M = 2 Chern Simons matter theories 

In [19] it was shown, by using localization techniques first introduced in this context in [35], 
that the partition function of N > 3 Chern-Simons-matter theories can be written as a matrix 
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integral. This result was generalized to = 2 theories in [14, 9]. We now give a brief summary 
of the ingredients involved in these N = 2 matrix models. A basic building block is the function 



l{z) = -zlog (1 - e 2 ™) + - (ttz 2 + -Li 2 (e 2 ™) - 



1 



7T 



17T 

12' 



(2.1) 



which has the following basic properties: 
1. It is odd 



£(z) = -£(-z). 



(2.2) 



2. Using the standard expansions of the logarithm and dilogarithm, one immediately shows 
that £{z) has the asymptotic expansion 



£{z) = ± 



17T 



m=l 



m 2txw? 



(2.3) 



where the ± sign corresponds to Im(z) > 1 or Im(z) <C — 1, respectively. 
3. It satisfies the equation 



d£ 
d~z 



-TTZ COt(7Tz). 



4. If we denote 
and 

we have from (2.3) 



z = A + 



2tt 

A = 1 — A 



e(z)+£(z*) = -A\e\-v A (6), 
where the function va(0) has the following expansion 



VA 



-m|0| 



m=l 



(R 



;j | sin ^27rmA^ cos ^27rmA^ 



Trm Txm^ J \ ) m \ 

In this paper, this function will be interpreted as a short-range potential. 
5. With the above notations, if A = 1/2, we have 

01 



(2.4) 

(2.5) 

(2.6) 
(2.7) 

|0| > 1. (2.8) 



£{z) + £{z*) = -log 



2 cosh ■ 



(2.9) 



6. It is not difficult to show that the Fourier transform of va(x) is given by the simple ex- 
pression, 

/27T 



55a M = 



1 - A sinh(2vr(A - l)w) 



IjJ TTUJ COsh(27Tw) — 1 

where the Fourier transform of a function f(x) is defined as 



1 f°° 
/(«) = -= / /(x)e 1 " 



c dx. 



(2.10) 



(2.11) 
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7. Let f(x) be a function. Up to exponentially suppressed terms, we have 

/oo / poo \ 

dx'v A (x - x')f(x') » I / dx'^A^') /(x). 
-oo \J — oo / 

The integral appearing above can be computed by using the Fourier transform (2.10): 

2tt 2 



(2.12) 



f 



dx va(x) 



2t™ a (0) = ^-(A - 1) (2A 2 - 4A + 1) . 



(2.13) 



We now consider a general N = 2 quiver CSM theory, made up of nodes connected by edges. 
Each node has a U(N) Chern-Simons gauge theory with level k a , where a = 1, ■ ■ ■ , r labels the 



nodes. Each node has an associated to it set of N eigenvalues A. 



(a) 
i ' 



1, • • ■ ,N, and the matrix 



integral is obtained by integrating over all the eigenvalues. We will now list the contribution 
of the different fields to the integrand of the matrix integral. We denote by Aj the variables 
corresponding to the a node, and by ^ those corresponding to the b node. 
The contribution of the CS vector multiplet at the node a gives a factor 



n(2sinh^ 



i<j 



(2.14) 



while the classical CS action leads to 



iA:, 



N 



cxp 



47T 



(2.15) 



i=l 



We will assume, as in [15], that the nodes are connected by pairs of bifundamental fields A a i>, 
Bb a , with anomalous dimensions A( a ,b) an d A(b,a)i respectively. This leads to the following factor 
in the integrand 



A r 



Ai - n 3 \ (~ \j - nj 



(2.16) 



A field in the adjoint representation in the a-th node is represented by 

, Aj — Xj 



JJexp 

ij 

while a field in the (anti) fundamental gives 



M A a + i- 



2vr 



(2.17) 



n exp 



2vr 



(2.18) 



2.2 The interacting Fermi gas picture 

In [28] it was shown that, in the case of = 3 Chern-Simons-matter theories, the matrix integral 
obtained from the localization approach of [19, 14, 9] can be re-expressed as the partition function 
of a quantum, one-dimensional, non interacting Fermi gas. The number of fermions, N, is simply 
the rank of the gauge group U(N). Let us review this result by using a formalism suitable for 
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generalizations. We will consider the generalization of ABJM theory given by necklace quivers 
with r nodes[12, 18], and with fundamental matter in each node. These theories have a gauge 
group 

U(N) kl x U(N) k2 x • • • U(N) kr (2.19) 

and each node will be labelled with the letter a = I,--- , r. There are bifundamental chiral 
superfields A aa+ ±, B aa -i connecting adjacent nodes, and in addition we will suppose that there 
are Nf a matter superfields (Q a ,Qa) in each node, in the fundamental representation. We will 
write 

k a = n a k, (2.20) 

and we will assume that 

r 

J> a = 0. (2.21) 

a=l 

According to the rules reviewed above, the matrix model computing the S 3 partition function 
of this necklace quiver is given by 



2 sinh 



in a k \ 2 r 1 I 



W J - ^ (2cosh^)^i i i n„2cosh(^±M) " 



a.; 



The building block of the integrand in (2.22) is the following A-dimensional kernel, associated 
to an edge connecting the nodes a and b: 

N TT 2 sinh ( A '~ Aj ^ 2 sinh [ ^ 

K ah (X^ ...,X N ;^,..., m ) = -Lj[ e -v-M nt<3 - \ \ J ,^ * (2.23) 



i=i 



n u 2coflh(*=fii) 



Here, 



U a (\) = -^A 2 + N fa log ( 2 cosh (2.24) 

and it will be interpreted as a one-body potential for a Fermi gas with N particles. 

To make the connection with a Fermi gas, we introduce the projection operator on totally 
antisymmetric states 

P= M £ (-i) e(ff) "> ( 2 - 25 ) 

which satisfies 

P 2 = P. (2.26) 

Let 

|Ai,---,Ajv) (2.27) 

be the basis of space eigenstates for an A-particle system Hn- We introduce the appropriately 
antisymmetrized states 

|Ai, • • • , Aat} = P|Ai, ■ ■ ■ , \ N ) = ^ (-^l^d)' ■ ■ ■ > X «(N)) (2-28) 
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which are a basis of the Hilbert space of fermions J 7 ^ (see chapter 1 of [32] for a very useful 
summary of these properties). We now want to interpret the kernel (2.23) as a matrix element 



K a b(M, ■ ■ ■ , Atv; hi, ■ ■ ■ ,hn) = {Ai, • • • , Xn\ Pab • • • ,Pn} , 



(2.29) 



in terms of a non-symmetrized density matrix p a t> (i.e. a density matrix for distinguishable 
particles). We first notice that 

{Ai, ■ ■ ■ , Ajv| p ab ■ ■ ■ , p N } = "^y ^2 ( _1 ) £(C7) Pa6 (Ai, ■ ■ ■ , X N ; ^(i), • • • , Ma(AT)) • (2-30) 



We now use the Cauchy identity 



n 



i<j 



2 sinh 



2 



2 sinh 



n M 2cosh(^) 



det 



2 cosh (^) 



(2.31) 



2 cosh 



In this equation, Sn is the permutation group of N elements, and e(a) is the signature of the 
permutation a. We obtain, 



K ab (\i,--- ,X N ;pi,--- ,p N ) = TTe ^Wty ( -\ 

Nl f = \ \2cosh^2^/ 

n iV AT 



(2.32) 



where we denoted 



t(x) 



2 cosh | 



By comparing with (2.30), it follows that 



N 



N 



Pa&(Ai,-" ,Aiv;/xi,--- ,/xjv) = JJe JJ *(Ai - ^). 



(2.33) 



(2.34) 



This factorization tells us that the iV-particle system is non-interacting, since the density matrix 
is completely factorized. 

A more general construction is possible which takes into account further interactions between 
the nodes. In fact, any kernel K which is antisymmetric in \, fj,j defines a density matrix p 
through the equation 



K (Ai, ■ ■ ■ , Xn; pi, ■ ■ ■ , Pn) = {Ai, • • • , \n\ p \p\, ■ ■ ■ , pn} 
Let us then assume that we have a more general kernel, of the form 
K a b(Xi, ■ ■ ■ , Xn; pi, ■ ■ ■ , Pn) 

1 



(2.35) 



-/ab (Aj,Mj) 



(2.36) 



h3 



2 cosh 
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which includes a two-body interaction in the node b, given by the potential Wj,, as well as a 
general interaction between the two nodes a, 6, f a t,. We first notice that 



n- 



n- 



(2.37) 



are invariant under any permutation of the Aj, p-i (since we are taking a product over all possible 
pairs of particles). Therefore, the r.h.s. of (2.36) defines an iV-dimensional kernel p ab . We can 
write 



N 



Po6(Ai,--- .Ajv^i,-" ,^) = ]le- u '^]le-f^^]le- w ^^l[t(X i -^). (2.38) 



i=l ij ij 

We will assume, as it happens in the N = 2 matrix integrals, that 

/afe(A,/i) = fab(\ - /i), W 6 (//,//) = W 6 (// - fjf). 

We can then write, 



(2.39) 



N 



Pab(M,- m - ,Ajv;Mi,-- - > Mat) 



n e_C7 ° (A<) 5a 6 (Ai - /*) JJe-^^-w) JJe-^^-w), (2.40) 



i=i 



where 



9ab(y) = e-f<*<*h(y) 



(2.41) 



incorporates the diagonal elements i = j in /„;,. 

Let us now suppose that we have an N = 2 matrix integral associated to a quiver with r 
nodes, where the CS levels are given again by (2.20). We use the Cauchy identity (2.31) for each 
node of the quiver. It is then easy to see that the one-body potential associated to the node is 
given by 



U a (X) 



in a k 
Att 



X 2 -N f J[ A fn +i 



"2tt 



NjJ f A 



. A 

i; - 1 

f* 27T 



(2.42) 



where Nf a (Nj ) is the number of (anti) fundamentals in the o-th node. If there are fields in the 
bifundamental connecting the edges a and b, the function f a b is given by 

f ab (y) = -£ (A (0)6) (A (6)0) - i^) - log (2 cosh I) , (2.43) 

so that the function introduced in (2.41) is 



g ab (y) = exp [l (a (c1)6) + + * (a (M - 



(2.44) 



Fields in the adjoint representation contribute to the function that we have denoted by Wb(ni — 
fij), which is then given by 



Y f w h { JH -nj) = m (A a )+XJ 



e A a + i 



. Hi - Hj 



2vr 



+ i Aa-i 



. Mi - Mj 



2tt 



(2.45) 



As in [28], we identify p a & as a canonical density matrix, defining a Hamiltonian for a Fermi 
gas. To write this Hamiltonian in a more explicit form we use the Wigner transform. We recall 
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that the Wigner transform of an operator A is given by (see [11, 39] for a detailed exposition of 
phase-space quantization) 



A 



(2.46) 



The Wigner transform of a product is given by the ^-product of their Wigner transforms, 

(Ab) w = A W *B W , (2.47) 



where the star operator is given by 



★ = exp 



~2 



and 



h = 2irk. 



(2.48) 



(2.49) 



It is convenient at this point to rescale the variables as A — > X/k. The Wigner transform of the 
canonical density matrix is then given by, 



PZ = e" ^ U ^ /k) * F ab ({q t }, {p.}) * e" E « , 



(2.50) 



where 



dyi • ■ ■ dy N e^=i ~s~ _[ J_ exp <j -/, 



Qi ~ Qj l _ Vi± Vj 
k 2 k 



N 



Ylgab(yi/k). 

(2.51) 



i=i 



This can be also written as 



N 



i=l 



^(te},ta})= / lldy; ,/,,,!//, /,•)(• J]' 



i=i 



^"«i; exp< j_ /a6 



1i - qj 



Qi ~ Qj 
k 



where 



/ipy 
dyg a b(y/k)e~*-. 

The above construction defines the iV-body Hamiltonian associated to p a b as 



Pab — e * 



(2.52) 

(2.53) 
(2.54) 



The total density matrix associated to the quiver is simply obtained by taking the product of 
the p a b over all edges. This is based on the identiy 



dA« • • • dAMtfx (A«, • • • , A«; xf\ • • • , A«) ■ ■ ■ tf r (a<'\ ■ ■ ■ , A«; A?\ • • • , A«) 



• • A {1) 



PKi ■ ■ ■ PX r 



^i V ■ ■ '^i^}' 



(2.55) 



- 8 - 



which just follows from using the resolution of the identity in Fn, 



dX | Ai, • • • ,Ai}{Ai,--- , Ajv| = 1, (2.56) 



r — 1 times. The corresponding, total Hamiltonian is simply obtained by taking the star product 
of the Wigner transforms p^ b . For example, if we have a circular quiver, we obtain 

e/ w =e, w * e* w * ■ ■ • * e* w * e* w . (2.57) 

In theories with N > 3 supersymmetry and canonical anomalous dimensions A = 1/2, 
the resulting Fermi gas is a non-interacting one, as shown in [28]. We then obtain a one-body 
Hamiltonian whose kinetic term, at leading order in h, is given by 

- log g ab (p) = log (2 cosh ^) . (2.58) 

For a general quiver, the Hamiltonian one obtains from the above procedure is quite complicated. 
At leading order in h, we have a one-body Hamiltonian whose potential term is given by (2.42), 
summed over all nodes, and a kinetic term given by — log g a b(p), summed over all edges. The 
functions / a &, Wt, lead to two-body interactions. The function Wj, leads to a standard interaction 
between the fermions, while the function f a b leads to a more complicated, non-standard inter- 
action between fermions which is velocity-dependent (since its Wigner transform involves both 
the positions qi and the momenta p.j). In addition, when considering h corrections, one will get 
fc-body interactions from the commutators of these one-body and two-body interactions. As we 
will see in the next section, an important simplification occurs in the theories with one single 
node. Before looking at this case, we will analyze the fermion system at large N and see how 
the functional theory developed in [10, 15] appears as a mean-field theory for this Fermi system. 

2.3 Large N limit and Thomas Fermi approximation 

The calculation of physical quantities in an interacting Fermi system and in the presence of an 
external potential is a non-trivial problem, and one is usually led to approximation schemes. One 
such scheme, which in some cases becomes exact where the number of particles ./V is very large, 
is the Hartree approximation and its semiclassical limit, the Thomas-Fermi approximation. Let 
us quickly review some ingredients of this approximation, in the simple case in which we have a 
one-body Hamiltonian of the form 

h(q,p)=T(p) + U(q) (2.59) 

and a two-body interaction with potential V(q, q'). The Hartree approximation can be regarded 
as a saddle-point evaluation of the many-body path integral (see for example [32, 21]). This 
saddle is characterized by a mean-field density p(q) and an effective one-body potential 

U c $(q) = U(q) + J dq' V(q, q')p(q'). (2.60) 

The second term in the r.h.s. is the Hartree (or direct) term. The density p(q) is detemined in 
a self-consistent fashion by the Hartree equation (at finite temperature) 

P(q) =n heB (q,q), (2.61) 
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where 

is the matrix element of the average number operator, p is the chemical potential, and the 
effective Hamiltonian h e s is 

h cG (q,p) = T(p) + U cS (q). (2.63) 
The density satisfies the normalization condition 

dqp(q) = N(p). (2.64) 

which gives the relation between N, the number of particles, and the chemical potential p. From 
this one can compute the grand potential in the Hartree approximation, 



dJ 
dp 



dqp(q), (2.65) 



and therefore all the thermodynamic properties of the system. 

In general, the Hartree equation (2.61) is not easy to solve. However, one can do further 
approximations. First, we note that the diagonal matrix element can be evaluated in terms of 
Wigner transforms, 

dp w 



nH{q,q) = J ^^n H {q,p). (2.66) 
On the other hand, in the semiclassical limit, we can use the classical Hamiltonian, 

n fffaP>* eftH(J-M) + l ' (2 - 67) 

and in this limit (2.61) reads 

P(q) = J ^ [exp (0 (T{p) + U cS (q) - p)) + If 1 . (2.68) 

This is the Thomas-Fermi equation at finite temperature, see for example [37]. A further ap- 
proximation involves going to zero temperature. Then, the semiclassical occupation number is 
determined by 

n^(q,p)^Q{p-H(q,p)). (2.69) 

The "classical" Fermi surface, 

H(q,p)=p, (2.70) 

defines implicitly the so-called local Fermi momentum 1 pp(q,p). The density p(q), in the 
semiclassical limit and at zero temperature, is called the Thomas-Fermi density. In the zero- 
temperature Thomas-Fermi limit, (2.61) can be regarded as an extremization condition for a 
density functional, the Thomas-Fermi functional, 

E TF [p] = t TF (q) + J dqp(q)U{q) + l - J dqdq'p(q)V(q - q')p(q') - p {^J dq p(q) - N^j , (2.71) 

1 In this part we assume that the Hamiltonian is symmetric in momentum. 
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where £tf((z) is the kinetic energy functional, 



*tf(?) = J ^ T (P)0 " h cS (q,p)) . (2.72) 

Notice that, in one dimension, the Thomas-Fermi density p(q) is always proportional to the local 
Fermi momentum, 

P(q) = ^PF(q,n), (2.73) 

where in the r.h.s. we have taken the positive solution for the local Fermi momentum. 

Of course, for an "atom," i.e. a three-dimensional Fermi gas in an external, attractive 
Coulomb potential, and with mutual Coulomb repulsion, the above formalism leads to the stan- 
dard Thomas-Fermi approximation in atomic physics. This approximation become exact as 
A — > oo [23]. As an even simpler example of the Thomas-Fermi formulation, we can consider 
the non-interacting Fermi gas which appears in AB JM theory [28] . In AB JM theory, the Hamil- 
tonian is, at leading order in h, 

H(q,p)= log (2 cosh |) + log (2 cosh |) « M + M . (2 . 74 ) 
The Thomas-Fermi kinetic functional is 

tTFiq) = I ^h4 e{11 - H(q ' p)) = T p2(q) - (2 - 75) 

Therefore, the Thomas-Fermi functional is simply (up to the constant pN) 

Etf[p] = J dq 



(2.76) 



Extremizing this functional leads to 



Pfe) = (2.77) 



which is just (2.73) for the Fermi momentum obtained from the Hamiltonian. The constraint of 
having A particles gives 

(i = — - — (2.78) 
and the evaluation of this functional on the above solution reproduces the well-known result [6] , 

E TF [ p ] = ^Kk^N 3 / 2 . (2.79) 

The general Hartree/Thomas-Fermi approximation is supposed to be exact, at large N, if the 
system becomes very dense in that limit (see [38] for some useful comments on this issue). We 
see from (2.77) that the size of the Fermi gas scales like A 1 / 2 , therefore the density grows like 
A 1 / 2 at large A, and indeed we are dealing with a dense system. This behavior is in fact typical 
of all the systems considered in [28] and in this paper, so the Thomas-Fermi approximation is 
appropriate. 

The functional (2.76) is very much like the one in [10] after the y variable has been integrated 
out. The factor p 2 , which in [10] arises from the interaction between eigenvalues, is due here to 
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the kinetic term in the Thomas-Fermi functional. We will now show that a mean-field treatment 
of the general interacting Fermi gas considered above leads to a functional of the Thomas-Fermi 
density p which is identical to the one obtained in [15, 29, 4]. 

First, let us note that there is a freedom in the choice of p in (2.35) for a given K. In 
principle one can simply choose 

Pafe(Ai, • • • , Ajv; pi, ■ ■ ■ , hn) = K ab (Xi, • • • , A^; p±, ■ ■ ■ , pn)- (2.80) 

Let us separate in this kernel the one-loop and the classical contribution due to the CS action: 

k 2 

K ab ({\i}, {pi}) = e^ K a b({Xi}, {pi}). (2.81) 
The Wigner transform of the one-loop part is given by 

*S««>.ta» = /nw^*-({f + £}.{|-&})- P.82) 

i 

According to the general principles of the Thomas-Fermi approximation, one can suppose that 
the fermions are distributed uniformly with the density l/(27i7i) in the domain of the phase 
space bounded by the upper and lower Fermi momenta p~^{q) and Pp(q). Then the density of 
the distribution in the coordinate space is given by 

Let us note that in (2.82) y; L plays the role of the difference of the eigenvalues associated to 
adjacent nodes. Suppose that the main contribution is given by imaginary m (after deformation 
of the contour of the integration). Then, by taking a continuous limit, 

Vi -> iVabiQ), (2-84) 

the kernel K ab can be computed by repeating the calculations done in Appendix A of [15] (see 
also [17, 16, 10]). It has the following form (for the sake of simplicity we consider the case without 
fundamental matter multiplets): 

(If + II • If " ID M exp {/ d "" 2 (")^(«("))} • < 2 - 85 > 



where 



Also 



where 



My) = fc/A (ai6) (2^) + fc/A (6 , a) (-^) + kf Aa (0), (2.86) 
2vr 2 

Ia(z) = —(z + A)(z + A- l)(z + A -2), 0<A + z<2. (2.87) 

^Yl PiVi ~ ~2^fr? J ' dq / dppVab ^ = ~\ J ' d 1P(l)P(l)yo.b{Q), (2-88) 

m s ?M±i£M. (2 . 89) 
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The Wigner transform of the total density matrix is given by 



p W (fe},fe}) = i K a+1 ({Qi}dPi}) = i e^^*i^ w a+1 (fe},fe}) 

o=l a=l 



= ? K^ +l ({ qi },{ Pt -X)n 6 ffi}). (2.90) 
a=1 6=1 

Then 

/> W (fe},fe}) ~ J y dgp(g) I ® a<1 - P ^ ■ y a ,a+i(q) + p(g)0o,o+i(3/o,a+i(g)) j , 

(2-91) 

where Q a = X^6=i n 6 an( l we replaced the ^-product by the ordinary one. Using the saddle 
point approximation for the integration over y a , a +i we arrive at the following expression for the 
Thomas-Fermi functional: 



Etf[Pf,P f ] = min V f dqp(q) \ ® a<1 P ^ • y a ,a+l{q) ~ p{q)(t) a ,a+l{ya,a+l{q)) 

a v 

dqp(q)-N^j , (2.92) 

where p(q) and p(q) are related to the Fermi momenta p~^{q) andpp(g) through (2.83) and (2.89). 
The thermodynamic limit of the free energy is given by 



F= min E TF \pJ,,p F ] 



_ mm 



J Mp{q) | ^ ag • y ata+1 (q) - p(q)4>a,a+l(ya,a+l(q)) 



-H[ / dqp{q)-N 



. (2.93) 

This minimization prescription is the same as the one that was first introduced in [10]. The 
function p can be considered as the Lagrange multiplier for the condition J2 a Va 0+1 = (cf. 
[17])- 



3. Theories with one single node: flavored theories 
3.1 The Thomas Fermi approximation 

We consider now the flavored theories with one single node studied in [3, 14] and whose matrix 
model was analyzed in section 4 of [15]. These theories are much simpler from the point of view 
of the picture in terms of interacting fermions, since at leading order in the K expansion they 
have only position-dependent interactions. 

The theories under consideration have a single node with gauge group U(N), and three sets 
of pairs of chiral superfields in the fundamental representation, 

(qf\qf ] ), * = i,2,3, j = i,2,---,m. (3.1) 
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Their anomalous dimensions are denoted by 



There are also three adjoint chiral superfields 

X±, X2, X3. 

The anomalous dimensions of these fields are denoted by 

A,, i = 1,2,3. 

They satisfy the following constraints 

X> = 2, 

and 



i=i 



A 9l + A 5 . + Aj = 2. 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



One has to add monopole operators [15], but as shown in that paper they do not contribute to 
the final answer for the free energy and we will not consider them. Their inclusion is however 
straightforward. We will introduce a parameter k as follows, 



n i — kfi, 



(3.7) 



and we will formally regard it as the parameter for a semi-classical expansion in the Fermi gas, i.e. 
the Planck constant is given by the same relationship as (2.49). The matrix model computing 
the free energy on the three-sphere of this theory can be obtained with the rules reviewed in 
subsection 2.1. Its integrand contains the factors 



nn«p 

i=l m,n 

x JJexp 



1 " A; + i 



Am A n 

2^ 



. A m \ 



1 A (/ , -,/!l-A,-i^ 



2vr J 



2ir J 



(3.8) 



due to the adjoint superfields and the fundamental superfields. In addition, we have the standard 
contribution 



n 



2 sinh 



(3.9) 



due to the U(N) vector multiplet. We now introduce the quantity 



Q({M) 



U m <n (2 sinh 



2 J 



n m ,n 2c ° sh 



1 



(3.10) 



o-eSjv 1 1m ^ \ 2 

The interaction between eigenvalues in the one-node matrix model can then be written as 

a(Ar)Q ({|}) n &~ v{qm ~ qn \ (3.11) 
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where 

\m = ^, (3.12) 

the potential is given by 

V( q ) = - ± [« (2L, + i JL) + « (A, - i JL)] - 2 log (2 cosh 4) (3.13) 

1=1 

and the prefactor in (3.11) is 

a(iV) = exp|iV^log2 + ^^(A i )^|. (3.14) 

It comes from the diagonal terms m = n not included in (3.11). 

After using the Cauchy identity, the integrand of the matrix model can be written as 

{Ai, • ■ ■ , Ajv| p I Ai, • • • , Aat} , (3.15) 

where 

N N 

P (Ai,-- - ,\n;m,--- ,m) = n e ~ u(kXm) n e- y ( fcA ^ feA -) n ^ 

m=\ m<n m=l 

t(y) is given in (2.33), and the one-body potential can be read from (3.8) and it reads, 

U (q) - - ±n, [l (l - A ai 4- i JL) + / (l - A. - 1 JL)] . (3.17) 

1 = 1 

The Wigner transform of p is given by 

N N N 

i»N{{qi},{Pi}) = n e ~ hu(qm) n *-* nqm - qn) * n 2cosh ^ * n e " |c/(,?m) n e-^-^ 

m=l m<n m=l \ 2 / m=l m<n 

(3.18) 

Up to corrections coming from higher orders in the star product, this is a gas of N particles with 
kinetic term 

T(p)=log(2cosh^), (3.19) 
as in [28], and a one-particle potential U(q). At large \q\ this potential becomes 

tf(?) = 7y + ■■■ , (3-20) 

where 

3 

1 = Y J h^- (3-21) 

i=l 

The interaction potential between the particles in this gas is 

F W = X>, (f ) - 2, 1/2 (f ) , (3.22) 

i=l 
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since, due to (3.5), the long range potential proportional to \q\ cancels. The resulting potential is 
a repulsive, short-range potential. Of course, on top of the one-body and the two-body potential, 
the star product leads to an infinite series of "quantum" corrections involving s-body poten- 
tials, which are in addition velocity-dependent (they depend on both qi and pi). The quantum 
Hamiltonian Hy/ can be computed by using the Baker-Campbell-Hausdorff formula. One finds: 

#W(M, {Pi}) = Yl U ^n) + E T ^«) + \Y. ~ Irn) 

n n rriy^n 

~ T2 £ ( T/( ^ } ) 2 u " {qn) - T2 £ T ' {Pn) (T ' (Pn) " T ' {Vm)) v " {qn " qm) 

+ DE T '>™) \ U'{qn)+Y. V '^-q m )\ +0(h 4 ). (3.23) 



We will now analyze the resulting Fermi gas after doing the following approximations: 

1. We neglect the quantum corrections to the Hamiltonian. 

2. We treat the two-body interaction in the Hartree/ Thomas-Fermi approximation reviewed 
above. 

3. We take the limit of zero temperature. 

4. We take the polygonal limit of the one-particle Hamiltonian, corresponding to large \q\ and 

\p\- 

As we will see, we will recover in this way the solution in [15] by considering the so-called 
Hartree approximation to the interacting problem and doing moreover the following approxima- 
tions: 

In the polygonal limit, the kinetic term is simply given by 

Tip) « ^ (3-24) 
and the zero-temperature Thomas-Fermi equation reads in this case 



p(q) 



J 2% @ - l 4 - u{q) - 1 dq v{q - q)p(q) ) • (3 - 25) 



In one dimension, the Thomas-Fermi density is just given by the local Fermi momentum, 

M = (3-26) 

where PF{q) solves the equation 

PF {q) = 2 M - 217(g) - 2 J dq' V(q - q')p(q'). (3.27) 



- 16 - 




Figure 1: The Fermi surface for the interacting Fermi gas associated to flavored one-node theories, in 
the Thomas-Fermi approximation. 



In the polygonal limit, in which we neglect exponential corrections to the functions appearing 
here, we can use (2.12) to write 

J dq' V(q - q')p(q') « (J dq' V(q')\ p(q). (3.28) 

Using also (2.13) and (3.5) we find, 

j dq' V(q - q'W) « ir 2 k (5 - 1) p, (3.29) 

where 

5 = 4AiA 2 A 3 . (3.30) 
After all these approximations, the solution to the Thomas-Fermi equation is simply 

P(«) = ^L? (3-31) 
and we can use the polygonal approximation (3.20) for U(q). The effective potential is 

U eS {q)= (l-^\p+^U{q). (3.32) 

In the Hartree approximation, we can think about the interacting problem as a one-body 
problem with the effective potential above. The Fermi surface is then determined by 



or equivalently 



y + U cS (q) = fi, (3.33) 



2 + 5 2 " 5' [6 - S ' 
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and we have depicted it in Fig. 1. The area of this Fermi surface, measured in units of 2irh, 
equals the number of particles of the gas, 



N = 



vol(/i) 1 2/i 2 



2ttH 5j ir 2 k 



(3.35) 



This is of course equivalent to the normalization condition (2.64). The grand potential is, at 
leading order, 

After an inverse Legendre transform, 

F(N) = J(ji(N)) - N»(N), 

we immediately find 



(3.37) 



F(N) « -^LN^kWy/Sj = - 2 -^N^ 



\ 



A ± A 2 A 3 [^niAi 



(3.38) 



,i=l 



which is exactly the result of [15]. 

Notice that, in the mean-field picture coming from the Hartree/Thomas-Fermi approxima- 
tion, we have iV fermions in the presence of an external, linear confining potential, and they fill 
out an interval of lenght ~ fi ~ iV 1 / 2 , just like in ABJM theory. The effect of the short-range 
interaction is to modify the precise numerical value of the parameters involved in the solution, 
without changing the qualitative picture of the Fermi droplet. Since we have a very dense sys- 
tem at large N, we expect the Thomas-Fermi approximation to be exact in the large N limit. 
Moreover, since the support of the Thomas-Fermi density grows like iV 1 / 2 , we are justified in 
taking the polygonal approximation, since corrections to this approximation will be subleading, 
as in [28] (as we will see in a moment, in the presence of a long-range attractive potential, the 
support of the Thomas-Fermi density will be of order 0(1), and the polygonal approximation 
breaks down.) 

We can then interpret the functional obtained in [15] as the Thomas-Fermi approximation 
to an interacting Fermi gas, leading to the correct result in the large N limit. In particular, 
(3.36) is the correct grand potential at leading order in fi. 

3.2 Corrections to the Thomas Fermi approximation 

One advantage of the Fermi picture, already emphasized in [28], is that one has in principle a 
systematic way of improving the leading large N approximation. In the non-interacting case 
considered in [28], this made possible to determine an infinite number of subleading 1/N correc- 
tions. The interacting case is much more difficult, but we expect the following corrections to the 
above result for the grand potential (3.36): 

1. As we already pointed out, without leaving the Hartree approximation, there will be cor- 
rections coming from finite T effects and from deviations from the polygonal limit. These 
corrections already occur in the context of the ideal Fermi gas analyzed in [28] , and there 
are clearly present in this example as well. In order to calculate these corrections quanti- 
tatively we have to solve the Thomas-Fermi equation for p{q) beyond the polygonal limit 
approximation. 
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Figure 2: The exchange correction due to the two-body interaction (3.22) (left) and the first ring diagram 
contributing to correlation effects (right). 

2. One has to take into account the quantum corrections to the Hamiltonian. There are 
corrections to the one-body Hamiltonian, as in [28] , as well as to the two-body Hamiltonian. 
In addition, the quantum corrections will lead to s-body interactions, for all s > 3, which 
should be also taken into account. Furthermore, since we are using the Thomas-Fermi 
approximation, one should also expect corrections to the semiclassical approximation, of 
the Wigner-Kirkwood type. 

3. As it is well-known in many-body theory, the Hartree approximation is the starting point 
for a (resummed) diagrammatic expansion. Even if we restrict ourselves to the two-body 
interaction (3.22), there will be corrections coming from exchange and correlation effects. 

We expect on general grounds that these corrections will lead to sub-leading terms in \i 
in the grand potential (3.36). We don't have a systematic argument for this, but a detailed 
examination of many correction terms suggests that this is the case. More precisely, the next- 
to-leading correction to (3.36) comes for the exchange correction to Hartree theory, and goes 
like /j, 2 . However, the resulting correction to the partition function cancels against the prefactor 
(3.14). Let us see this in detail. 

The exchange correction to the Hartree approximation can be represented diagramatically 
as in Fig. 2, left (see for example [32]). Its contribution to the grand potential is given by 

JM = \j dqdq'n hcS (q',q)V(q-q')n hcS (q,q'). (3.39) 

We will evaluate this correction in the same approximation scheme that we used before, i.e. we 
will use a semiclassical calculation at zero temperature, and in the polygonal limit (improving 
this approximation leads to even more subleading corrections). The two-point function of the 
occupation number operator appearing in this expression can be computed in terms of Wigner 
transforms as 

n hcH (q',q) = J ^e^n&fop), (3.40) 

where 

R=^Y~> r = q-q'. (3.41) 
In the semiclassical and zero temperature approximation we have 

/ / x I dp *pr . , , 1 sin (pF(R)r/h) . , . 

n HeB (<l , 9 ~ / ^ eir0 ~ h cs(P, R)) = , (3-42) 

where pf{q) is the local Fermi momentum for the effective Hamiltonian. Therefore, the exchange 
term reads 

j M « JL f dRdr f ^MR)r/n) \\ { ^ (343) 
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Let us estimate this for large [/,. In the polygonal approximation, pp is given by 

PF(q) = f-^\q\ (3.44) 
and it scales as fi. We introduce now the rescaled coordinates £, ( as 

r = e//*, R = K (3-45) 

so that the exchange integral reads 

•Wrt - £ / dd? ( Si ^W )V(W, (3.40) 



where 



p F (R)=W F (0, MO = l-}\<\- (3-47) 
For large, we have then 

JWM « £v<o) / dCdj ( " fe | c|{/ " 1 )' = ^v(0) / dtf„(0 = iv(o W (3.48) 

where we used the normalization condition (2.64). However, we know from (3.13) that 

3 

F(0) = -2^(a/) -2 log 2, (3.49) 
i=l 

therefore the leading, fi 2 contribution of the exchange term cancels against the prefactor (3.14) 
in (3.11) (a similar cancelation occurs in mean-field many body theory, see section 7.2 of [32]). 

We have not found other sources for /x 2 corrections. It is easy to see that there are corrections 
of order 0(/x) (these corrections are already present in the non- interacting case considered in [28]) 
and of order 0(1). For example, the first ring diagram showed in Fig. 2 (right) can be evaluated 
in the semiclassical limit with the techniques of [21], and it can be seen to be of order 0(1). We 
then expect the grand potential to be of the form 

1 2// 3 

J ^ ) = ^zk +B » +A+ °^' (3 - 50) 

where B, A do not depend on fi. By using the standard inversion formula 

Z(N) = 7^-Jdfi exp [J(/i) - »N] , (3.51) 

much exploited in [28], we conclude that the partition function of the single- node theories is of 
the form 

Z(N) = C~ 1/3 e A Ai \c- 1/3 (N - B)} , (3.52) 
up to non-perturbative corrections in N. Here, 

This would confirm the conjecture made in [28] for this class of N = 2 theories. It would 
be important to perform a systematic analysis of the possible corrections to J(fJ-) in order to 
verify our preliminary analysis here. It would be also interesting to see if one can calculate the 
coefficient B by using the interacting Fermi gas picture. In principle, it is clear what one should 
do: among other things, one has to analyze the Hartree equation beyond the zero temperature 
and the polygonal approximation, and one has to take into account the corrections coming from 
many-body diagrams. 



- 20 - 



4. Theories with one single node: long-range forces 



We will now consider a different class of theories with one single node, namely N = 2, U(N) 
Chern-Simons theory with g adjoint multiplets. These theories were first studied in [8], and the 
matrix model partition function was analyzed in [30, 2] (for general g) and in [33, 31] (for g = 1). 
Here, we will study these theories from the point of view of interacting fermions developed in this 
paper. This will allow us to re-derive and improve some of the results in [30]. In particular, we 
will see how the Thomas-Fermi equation of this model - a single integral equation - determines 
the i?-charge of the adjoint multiplets in the large N limit and at infinite 't Hooft coupling in 
a very efficient way (the method used in [30] was based on the numerical extrapolation of the 
saddle-point equations at finite N and 't Hooft coupling). We also compute the large N free 
energy of these models from the Thomas-Fermi equation. 

The partition function on the three-sphere of j\f = 2, U(N) Chern-Simons theory with g 
adjoint multiplets is given by 



r N ■ f 

I \\&\ie^ Yl ( 2sinh 

i=l i<j ^ 



A?; — X', 



exp 



i,3 



1-h + i 



Aj — Aj 
2tt 



(4.1) 



This follows from the rules reviewed in section 2. We have denoted by h = A Ad the i?-charge of 
the adjoint hypermultiplet. In principle, this partition function depends on four parameters: N, 
k, g and h. However, the value of h is determined, as a function of N, k, and g, by maximizing 
\Z\ [14], so in fact there are only three free parameters. 

If we use the Fermi interpretation developed in section 2, we immediately find a long-range 
potential between the fermions of the form 



(g(l-h)-l)\x\ 



(4.2) 



Clearly, the nature of the fermionic system will depend crucially on the sign of the coefficient. 
If this sign is positive, we have a long-range attraction, and if the sign is negative we have a 
long-range repulsion. We will analyze these two situations separately. In both cases, as we will 
see, the Thomas-Fermi distribution is supported on an interval whose length does not longer 
scale with iV 1 / 2 , as in the M > 2 theories considered above and in [6, 10, 15, 28]. In the case 
where there is long-range attraction, the support of the distribution has a length of order 0(1). 
This means that, even at large N, we cannot use the polygonal approximation that we used 
in the case of a short-range potential, since the terms we would neglect in this approximation 
are as important as the terms that we would keep. Moreover, the quantum corrections to the 
Hamiltonian are also important in this case. 

The easiest way to incorporate all these corrections in the Thomas-Fermi equation is to take 
the standard large N limit directly in the matrix integral (4.1), in which the Thomas-Fermi 
distribution p(x) is interpreted as a density of eigenvalues. Standard techniques lead to the 
integral equation 



/ 



dyp(y)K(x -y) = n 



where 



a 



ax 



-ik 



x G supp p, 



(4.3) 



(4.4) 



and the kernel K is given by 



K(x) = gV h {x) - 2 log 2 sinh ■ 



(4.5) 
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where 

V»(») «(l- k + i i)-/(l-k-li). (46) 

In what follows we suppose that a is real and positive to ensure that the eigenvalues are dis- 
tributed along the real axis. The free energy for real values of k can then be obtained by analytical 
continuation. 

At large x, the kernel behaves as 

K(x) = (g(l-h)-l)\x\+0(e-M), (4.7) 

which is the long-range potential (4.2). At small x, it behaves as 

K(x) = -2log\x\+gA + Y,(A2k-^) ^y, (4.8) 

where A 2k = vf fe) (0) and B 2k are Bernoulli numbers. As we mentioned before, the distribution 
of eigenvalues depends drastically on the sign of (#(1 — h) — 1). 

4.1 The case of long-range attraction 

Suppose first that (</(l — h) — 1)>0. In this case there is a long range attraction force between 
eigenvalues. Let us rescale the density as 

p(x) = fif(x), (4.9) 

so that the integral equation (4.3) reads 



/ 



ex 2 



dyf(y)K(x-y) = l-—, (4.10) 



where 

e=— (4.11) 
and it vanishes at large (i. Then, from (2.65) we have 



N(ji) = dJ ^- = li J f(x)dx ee C M . (4.12) 



Therefore, in this theory, fi is proportional to N. This is again in contrast to the N = 3 theories 
considered in [6, 28, 10], as well as to the J\f = 2 theories with no long-range forces analyzed 
above and considered previously in [15]. Notice that the parameter e can be written, at large N, 
as 

where 

N 

A = - (4.14) 

is the 't Hooft parameter of the theory. The M-theory limit that we are considering here (N 
large, k fixed) corresponds as usual to the strongly coupled region of the 't Hooft parameter. 
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In particular, the limit of infinitely strong coupling A — > oo corresponds simply to e = 0. The 
canonical free energy in the large N limit is given by 

% = "MAO, (4-15) 
where the function n{N) is the inverse to N(fi), therefore 

f ^-^c n2 - ( 4 - 16 ) 

Qualitatively, the behavior of the interacting Fermi gas underlying this model is very different 
from the N = 2 Fermi gases with no long-range forces, as the flavored theory considered above. In 
the model with long-range attraction, the fermions, under the action of the attractive potential, 
form some sort of bound state whose size is of order 0(1) at large N. The system is still 
very dense, since the density grows linearly with iV, and therefore we expect the Thomas-Fermi 
approximation to give the right large N behavior. This is very similar to the analysis of baryons 
at large N in [38], whose size is of order 0(1) at large N but whose mass grows as 0{N). In 
our case, the Thomas-Fermi equation is simply the standard integral equation following from 
the large N density of eigenvalues of the matrix model (4.1). Let us note that the bound state is 
formed due to the long-range attraction, even if there is no external potential and e = 0. In this 
limit the integral equation has translation invariance which can be removed by assuming that 
the solution is centered at 0. 

The maximization principle of [14] says that, in order to compute h, we have to maximize 
— F, i.e. we have to minimize C. This determines the R-charge h of the adjoint multiplets. In 
the 't Hooft expansion, this R-charge, in the planar limit, is a non-trivial function of the 't Hooft 
parameter A. In our formalism, it is very easy to determine h in the planar, strongly coupled 
limit A — > oo (or e = 0), since the coefficient C\ e =o can be calculated numerically with high 
precision for any given g and h. To do this one first solves numerically the following integral 
equation 

' f g>h (y)K(x-y)dy = l, (4.17) 



/ 



obtained from (4.10) by setting e = 0. Then 

C\ e =0 = I fg,h{x)6x. (4-18) 



The large N, strongly coupled limit of h as a function of g was determined numerically in 
[30] by solving the discrete saddle-point equations for the matrix integral (4.1) at finite N, k, and 
then extrapolating the result to large N and A. The advantage of the integral equation (4.17) is 
that it gives directly the right R-charge h in the limiting region N, A — > oo without the need to 
do an extrapolation. 

An example of a numerical solution of (4.17) for the profile f 9: h is shown in Fig. 3. The 
behavior of C| e =o as a function of h for g = 2, 3 is shown in Fig. 4. Numerical calculations show 
that, for any real g > 1, there is always a minimum of C\ t= o at a certain h m i n (g) in the interval 
< h < 1 — 1/g defined by the condition that the long range forces are attractive. Therefore 

lim h min (g) = (4.19) 

which agrees with the result of [34, 30] saying that for g = 1 the R-charge of the adjoint field 
tends to zero in the strong coupling limit. Accurate numerical calculations show that /i m in(2) = 
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Figure 3: The graph of f g ,h(x), the solution of (4.17), for g = 2 and h = 0.2726. 
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Figure 4: The graph of C\ e= o as a function of h for g = 2 (left) and g = 3 (right). 



0.2726 ± 0.0001 and /i min (3) = 0.3539 ± 0.0001 (see Fig. 4) which is in agreement with the 
numerical values found in [30]. 

It was pointed out in [30] that analytic expressions for h could be obtained in the limit 
g — > oo. We will now show how to obtain a systematic expansion for h in this regime, from our 
integral equation (4.3). In this way we will recover and extend some of the results of [30]. Let 
us first make the following rescaling: 



£ C 
x = —=, y = — 

V9 y/9 



fix) = V fn ^ 



n>0 



Then (4.10) reads 



>ra>0 



^_ 2io g | { -;| + iogg + / 2 _^K-cg + 

9 9 V 9 J 2c? 



Let now 



(4.20) 



1-^-- (4-21) 
2ff 



(4.22) 



i=0 
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be the solution at order n, and 

C (n) = J / ( „)(0d£ (4.23) 

Then from (4.21) at £ = we have 

AoC\ n) = 1-/ dC/ ( ^„(0 (-^ + T + (^ - f ) £ + • • ■) ™ d f""" < 4 ' 24 > 
and after differentiating (4.21) w.r.t. £ we get 

'/^=f^^/«woO.-) («»> 

This equation can be solved iteratively: suppose we know the solution at order n — 1, /( n -i)(0- 
Then one can determine C( n ) from (4.24), and use (4.25) to solve for /(„) Let us note that for 
any given order n the r.h.s. of (4.25) is a polynomial in £, and thus (4.25) can be considered as 
the standard matrix model equation for fr n \ with a polynomial potential which can be solved 
explicitly and unambiguously once C7 n ) is given. 

To see how this works, let us compute the first few orders. At order we have the two 
equations, 

^4oC(o) = 1) 

[ f(0)(QdC _ e + A 2 /A ( 4 - 26 ) 
J C-C 2 5 

and the solution is 



/(0)(?H i±^_J__ e (4 , 7) 



At order 1 we have 



4>C (1) = 1 - - j dC/ ( o)K) (-2 log |C| + logs + A 2 Yj = 



1 1 



1 Ao 

log 5 + l + log(A 2 + ^ e) + 



A g [ &v z u ' 2A 2 + 4>e 

which can be rewritten as 



(4.28) 



l = A)g + logg + l + \og{A 2 + A Q e) + \ - ^ +o(-\ ■ (4.29) 
Using the procedure described above up to order 2, and specializing to e = 0, one finds 



1 

C 



- H<0), + .o g9+ | + log tf(0) - ^ (l - + O (i) . (4.30) 



£=0 

Minimizing C| e= o w.r.t. /i gives 



1 _ A _ 32 7 9 ) , „ ( 1 

2 7r 2 g 37r 4 <7" 



Vin = n ~ -5- - 7 + 73 ) • ( 4 - 31 ) 



The first correction to the asymptotic value 1/2 was already found in [30] with a related technique. 
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4.2 The case of long-range repulsion 

When 

-7 = (g(l -h)-l)<0 (4.32) 

there is a long range repulsion between the fermions (or the eigenvalues of the matrix model), 
and the external potential cannot be neglected at large N. After rescaling (as we will see later, 
/j, < in this regime) 

x = ^f ^ J^i y = V^V, p(x) = /(£) (4.33) 

the equation (4.3) becomes 

7 / f(ri)\t - V\dv + O (-^=) = 1 + ^. (4.34) 

Suppose /o(0 is the solution of this equation at large [i. After differentiating two times w.r.t. £, 
the equation (4.34) becomes 

27/o(0 = ( 4 -35) 

Therefore at leading order the distribution is constant in some interval (— £o 5 £o)> and zero outside. 
The end points can be determined by evaluating (4.34) at £ = 0: 



Then, 



47T7 1 




7 / -;— Md»7 = l => £o = \/ — ■ (4.36) 



and 



(4.38) 



F « — — 
3a 

Qualitatively, in the case of repulsion, the eigenvalues spread out at large N over an interval 
whose length grows linearly with N, and with constant density. Such configuration can be 
easily understood in terms of the forces between eigenvalues. Suppose we have some symmetric 
distribution, and let us consider an eigenvalue sitting at x. Let n(x) be the number of eigenvalues 
in the interval (— x, x). The long-range force between every two eigenvalues is constant and equal 
to 7. The total force acting on the eigenvalue at x from other eigenvalues equals 771(2;) and is 
directed outside of the interval (— x, x). This force should be compensated by the force from the 
external potential, which equals 

S ,4, 9 , 
and is directed towards the interior the interval (— x, x). Then 

jn(x) = ^ (4.40) 

and 

p(x) = = - — . (4.41) 

KW 2 dx 4tt7 v ' 
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The derivative of the free energy can be easily calculated as the energy of a probe eigenvalue 
added at the boundary of the distribution, and one finds 

The result (4.38) can be tested in the case g = 0. In this case, the matrix model (4.1) reduces 
to the CS matrix model of [24], and its free energy equals the free energy of CS theory on S 3 
with framing 1. The corresponding planar free energy reads (see for example eq. (4.31) in [25], 
to which one has to add i 3 /12 due to framing), 

F (t) = ^-^ + C(3)-Li 3 (e-*). (4.43) 
Here, t = g s N is the 't Hooft parameter and 

9s = "f. (4.44) 

Then, at large N and fixed k, 

which coincides with (4.38) in the case g = (7 = 1). 



5. Conclusions and open problems 

In this paper we have extended some of the results of [28] to the matrix models of N = 2 CSM 
theories, which we have formulated in terms of interacting fermions. The resulting system can 
be analyzed, at large N, in the Hartree/Thomas-Fermi approximation, and this leads to the 
formulation in terms of density functionals put forward in [10, 15]. Going beyond the large N 
approximation is in general difficult, although in the case of flavored theories with one node one 
can compute the next-to-leading correction by using the exchange correction to the Thomas- 
Fermi approximation. In the case of theories with long-range interactions, the Thomas-Fermi 
approximation gives an efficient way of calculating numerically the R-charge of adjoint multiplets. 

It is clear that this work leaves many open problems. One should find a more rigorous 
argument showing that the Thomas-Fermi approximation gives the leading contribution to the 
grand potential, and one should achieve a more detailed understanding of what kind of corrections 
are expected for J(/x), as it was done in [28] for the M = 3 theories. This might lead to a way of 
determining, at least in the flavored theories with one node, subleading terms in /j, in the grand 
potential, which would lead to a calculation of 1/N corrections to the anomalous dimensions. 
In the theory with long range forces, the Thomas-Fermi equation is nothing but the standard 
equation for the eigenvalue density, and one might try to go back to the traditional matrix model 
technology in order to determine 1/N corrections to the large N result. 
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